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We propose and demonstrate the method of population transfer by piecewise adiabatic passage 
between two quantum states. Coherent excitation of a two-level system with a train of ultrashort 
laser pulses is shown to reproduce the effect of an adiabatic passage, conventionally achieved with 
a single frequency-chirped pulse. By properly adjusting the amplitudes and phases of the pulses 
in the excitation pulse train, we achieve complete and robust population transfer to the target 
state. The effect is demonstrated experimentally by observing piecewise excitation of Rubidium 
atoms from 5si/2 to 5pi/2 electronic state. We show that similarly to the conventional adiabatic 
passage, the piecewise process is insensitive to the total excitation energy as long as the adiabaticity 
conditions are satisfied. The piecewise nature of the process suggests that robust and selective 
population transfer could be implemented in a variety of complex quantum systems beyond the 
two- level approximation. 
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The existence of robust and selective methods of ex- 
ecuting population transfer between quantum states is 
essential for a variety of fields, such as, precision spec- 
troscopy and atomic clocks [l], quantum computing [2], 
control of molecular dynamics and chemical reactions [3^. 
Traditionally, population transfer between two quantum 
states has been achieved by either executing a half-cycle 
Rabi flopping, i.e. by applying a "7r-pulse", or by op- 
tically inducing an adiabatic passage (AP) between the 
states of interest [5]. Though the application of a tt- 
pulse can be implemented on a very short time scale, it 
is far from being robust, as it is highly sensitive to fluc- 
tuations in the laser power, phase, and pulse duration. 
In contrast, AP, which for example can be executed by 
the STIRAP technique [6], or by "slowly" chirping the 
instantaneous frequency of the pulse 0, exhibits high de- 
gree of robustness against the fluctuations of many of 
the laser field parameters. Because of this property, AP 
with chirped pulses has been widely used for executing 
population transfers in atomic [8(P and molecular p^Of [TT] 
systems, including recent applications to enhancing the 
photo- association of ultracold molecules [T3] . 

Though well suited for two-level systems, population 
transfer with chirped pulses in multilevel systems be- 
comes sensitive to the exact value of the chirp and field 
strength [14^, T5^, thus losing some of its appeal as a ro- 
bust way of efficiently manipulating populations. When 
used with spectrally broad ultrashort pulses, pulse chirp- 
ing can no longer selectively populate a pre-chosen sin- 
gle state or a superposition of states JTJ |T5], and other 
methods of adiabatic [16 or non- adiabatic p!7| [TSj exci- 
tation must be employed. Since ultrashort pulses do 
not satisfy the adiabatic conditions, it is customary to 
treat the problem in a purely empirical way by design- 
ing feedback-controlled experiments with genetic search 



algorithms [191 HO]. However, these approaches lack the 
robustness and efficiency of the AP method. 

Recently, we have demonstrated theoretically that 
one can implement AP with ultra-short pulses by ex- 
ecuting the transfer of population in a piecewise man- 
ner[2l] |22l. The original work on Piecewise Adiabatic 
Passage (PAP) [21 proposed using two temporally over- 
lapping pulse trains. It was shown that by tuning the pa- 
rameters of the pulse trains one can bring about a com- 
plete population transfer between two quantum states 
through a third intermediate level. In this Letter we 
present the first experimental demonstration of the PAP 
method with a single pulse train by introducing the tech- 
nique of "piecewise chirping". We show that the piece- 
wise population transfer faithfully reproduces an AP pro- 
cess executed with a continuous frequency-chirped pulse, 
achieving the same level of robustness and efficiency. 

Consider a system of two states, |1) and |2), of energies 
El and £^2, respectively. An arbitrary coherent superpo- 
sition of these states can be written in the interaction 
representation in terms of two angles, and ^, as 

= cosO/2exp[-iEit] |1) smO/2exp[-iE2t ^ icj)] |2) . 

(1) 

This superposition state can be represented as a unit 
Bloch vector pointing along the (^, 0) direction. 

In conventional AP with chirped pulses, the system is 
driven by a near-resonant field with a time dependent 
frequency detuning A relative to the exact resonance fre- 
quency coq = {E2 — Ei)/h. The product of the field 
strength and the transition dipole matrix element can 
be conveniently expressed in terms of a Rabi frequency 
Q. Application of the Rotating Wave Approximation in 
a frame rotating at the transition frequency cjq, results in 
two dressed eigenstates, ip-^ and ip- . These dressed states 
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FIG. 1: (Color online). Two sample trajectories of the Bloch 
vector (thick gray arrow) during the AP process implemented 
with a single continuous chirped pulse (a), and a train of 
20 transform-limited pulses with quadratic temporal phase 
(b)(see text). 

constitute two stationary points on the Bloch sphere, 
(l9±,^±), where tan l9± = ±n/A, (f)± = J A{t)dt±7r/2. 
All other solutions of the time dependent Schrodinger 
equation precess around the axes defined by the above 
stationary points, displaying periodic trajectories corre- 
sponding to the familiar Rabi oscillations [5 . 

If the parameters of the field are varied adiabatically, 
and if initially A <C — 1^|, and at the end of the process 
A ^ there is an interchange of populations between 
state |1) and state |2). This interchange is expressed on 
the Bloch sphere as the movement of the Bloch vector 
from the 6 = direction to the = n direction and 
vice versa. The motion proceeds via the spiral trajectory 
shown in Fig. [ija). The above behavior holds as long as 
the adiabaticity condition is maintained [23^: 



\9m\ « V^Ht) + ^Ht) 



(2) 



We now consider the evolution of the Bloch vector un- 
der the action of a train of short, mutually coherent, laser 
pulses. Each pulse, centered at the resonant frequency coq 
and of duration r, generates a rotation P = R^(ap) of 
the Bloch vector about the y axis (7r/2,7r/2) by an angle 
ap = f^ft{t)dt. To account for the change in the carrier 
phase between two consecutive pulses we introduce an 
additional z-rotation on the sphere F = Rz{aF)' Thus, 
the overall evolution may be represented by a sequence 
of rotations U = ...FPFP... 

The product FP of two rotations of the Bloch vector is 
an overall rotation by an angle ao about an axis defined 
by the (6>o, ^o) angles, given to lowest-order expansion in 
ap, ap as 



ao 

00 

tan ^0 



= VI 



±7r/2-aF/2, 
iofp/ap. 



(3) 
(4) 
(5) 



By maintaining the same value of ap and ap through- 
out the pulse train we induce piecewise rotations of the 



Bloch vector around the closest stable point (6>o, 0o)- By 
slowly varying the values of ap,ap we can make the 
stable points move and the Bloch vector, captured near 
them, to follow faithfully. Intuitively, the conditions of 
such piecewise following are: (i) the y- and z-rotations 
should be small (i.e. that each pump pulse should induce 
an angular change much smaller than tt and each incre- 
ment in the carrier phase should be small too), and (ii) 
{^o^4'o) should move little during each pulse, i.e. 



AOo < V(ap2+ap2)/2. 



(6) 



If initially ap <C |ap|, the two stationary points, = 
and 00 = TT, correspond to the bare states |1) and |2). 
As ap and |ap| increase and decrease, respectively, the 
states originating in |1) and |2) follow a piecewise spi- 
ral trajectory towards the equator of the Bloch sphere. 
They cross the equator as soon as ap changes sign, and 
finally interchange with each other when ap decreases, 
approaching zero, while |ap| increases. An example of 
such a piecewise walk on the Bloch sphere is shown in 
FigQb). 

The above analysis shows that the adiabatic following, 
similar to that implemented with the continuous chirped 
pulse, can be executed by a sequence of pulses with slowly 
varying amplitudes, and with the absolute carrier phase 
changing from pulse to pulse in a non-linear way (i.e. de- 
creasing ap in the first half of the process and increasing 
ap in the second). An example of such AP, correspond- 
ing to the piecewise population transfer in atomic Kb, is 
shown in Figj2] Here, the field is given by a sequence of 
20 femtosecond pulses. For the k-th pulse {k G [0,20]) 



Ek{t) = Ak cos [uot + ^k] X sin^ 



t — tk 



(7) 



where tk marks the beginning of the /c-th pulse, with 
r = 300 fs being its full duration. The Ak amplitude 
represents the train envelope parametrized as a Gaus- 
sian of 3 ps width (fwhm). The frequency uoq is res- 
onant with the electronic transition between the states 
|1) = 5^1/2 and |2) = of Rb. The "piecewise chirp" 

of the pulse train is determined by the extra phase fac- 
tor ^k = — ^o)^/2, where ko = 11 and a = 0.2. 
As the pulse sequence proceeds, the pulse-to-pulse phase 
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FIG. 2: (Color online), (a) The amplitudes (oscillatory, blue) 
and phases (piecewise parabola, red) of the driving field, (b) 
The populations of states |1) (falling, blue) and |2) (rising, 
green) during the PAP process. 
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change, ap = OL{k — ko) + q;/2, smoothly evolves from 
large negative values to large positive values. Our simu- 
lations show that as long as the conditions of piecewise 
adiabaticity ([6| are maintained, the population transfer 
is robust with respect to the pulse shapes, their intensi- 
ties, and the exact value of the piecewise chirp a. 

The population transfer by piecewise chirping de- 
scribed here is a particular example of the Piecewise Adi- 
abatic Passage (PAP) concept introduced in Ref . [2T] . In 
the limit of infinitesimal rotations, the present process 
becomes equivalent to the usual continuous frequency 
chirping, associated with the quadratic phase change of 
the carrier phase in time. In this limit, the coarse-grained 
adiabaticity condition (|6| reduces, up to a numerical fac- 
tor, to the familiar quantum adiabaticity condition ([2|. 

In the proof-of-principle experiment, we studied the 
efficiency of the one-photon excitation of atomic rubid- 
ium (^^Rb) from the ground to the first excited electronic 
state, 5si/2 to 5pi/2, respectively. The effect of a single 
frequency-chirped laser pulse was compared with the ef- 
fect of a short pulse train described above. The pulses 
were produced by a standard Ti: Sapphire regenerative 
amplifier and had a spectral width of 7.1 nm (fwhm). 
We tuned the central wavelength of the laser to 795 nm, 
resonant with the 5si/2 transition. To generate 

a continuous or a piecewise chirp, the original pulse was 
spectrally shaped using a home made pulse shaper based 
on a double-mask liquid crystal spatial light modulator 
in the 4-/ configuration [24] . The excitation beam was 
focused with a long focal length lens (/ = 100 cm) onto 
a cloud of rubidium atoms continuously evaporated from 
the rubidium dispenser in a vacuum chamber equipped 
with the time-of-fiight ion detector. 

To determine the population of the excited state, we 
ionized the atoms by a second, probe pulse tuned to 1300 
nm and arriving in the chamber 2.5 ps after the excitation 
pulse. The ion signal provides good selectivity between 
the ground and excited state Rb atoms because the ion- 
ization of 5^1/2 requires two more probe photons than 
that of bpi/2' The power of the probe pulse was lowered 
to less then 0.1 fiJ so as to produce no detectable ions 
from the atoms in the ground state. To ensure ion sam- 
pling from the region of uniform excitation field strength, 
we confined the interaction region in the longitudinal di- 
rection by a 3 mm aperture placed between the rubidium 
dispenser and the laser beam (figjsj^a)). The suppression 
of the transverse spatial averaging was done by focusing 
the probe beam down to a smaller spot size than the size 
of the 795 nm beam (1/e^ beam diameters of 180 and 470 
/im, respectively). Finally, the effect of the laser power 
ffuctuations was eliminated by recording the energy of 
each excitation pulse together with the corresponding ion 
count. 

We used a home made spectral pulse shaper for in- 
troducing continuous frequency chirp and creating the 
pulse trains required for PAP. The resolution of the pulse 
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FIG. 3: (Color online), (a) Scheme of the experimental setup. 
Rb atoms first interact with a loosely focused pulse train (wide 
blue) and, after a short time delay (St = 2.5 ps), are ionized by 
a tightly focused probe pulse (narrow red). The total amount 
of ^^Rb^ is measured for each excitation pulse with a multi- 
channel plate detector (MCP); (b) Example of the applied 
spectral amplitude and phase shaping for generating a train of 
9 pulses with a piecewise chirp. Solid (blue) and dashed (red) 
lines show the amplitude and phase modulation, respectively. 
Dash-dotted line (black) represents the power spectrum of 
the original transform-limited pulse; (c) frog retrieval of the 
temporal field amplitude (solid blue) and phase (dashed red) 
for a pulse train of 9 pulses after applying the shaping shown 
in (b) Black dots correspond to the calculated quadratic phase 
of the pulses, All amplitudes are shown in arbitrary units. 

shaper of 0.17 nm per pixel allowed us to split the orig- 
inal pulse into a sequence of up to 9 well separated 
pulses. Continuous chirping was applied by means of the 
quadratic phase-only modulation, 0(a;) = a{u; — cjq)^, 
where a is the linear chirp and ujq is the center frequency 
of the pulse. Generation of a pulse train requires both the 
amplitude and phase modulation of the field spectrum. A 
sequence of N replicas of the original transform-limited 
pulse with the real amplitudes Ak and phases coot + 
separated by the time interval r were obtained with the 
following complex spectral mask: 

5K) = ^Afce-^(-"-+*'=)/|5|, (8) 

k 

where the frequency corresponds to the n-th pixel of 
the pulse shaper, and the normalization factor IS*! was 
used to satisfy \S{u;n)\'^ < 1 V n. Examples of the ap- 
plied spectral shaping and the corresponding pulse train, 
characterized by the method of frequency resolved optical 
gating (frog), are shown in Fig. [3]^b,c). 

In Fig. |4] we show the measured ion count correspond- 
ing to the population of the 5pi/2 state of rubidium for 
various parameters of the excitation field. The oscilla- 
tory (blue) curve in plot (a) represents the well known 
Rabi oscillations as a function of the energy of a single 
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FIG. 4: (Color online). Measured (dots) and calculated (lines) 
population of the excited state for a single pulse (a,c) and a 
train of 7 pulses (b,d). In (a) and (b), the efficiency of the 
population transfer is shown as a function of the integrated 
pulse energy. In both cases, the oscillatory data (blue dia- 
monds and solid lines) reflect the familiar Rabi oscillations, 
while the second set of data points (red circles and dashed 
lines) corresponds to the continuously (a) and piecewise (b) 
chirped excitation with a = 20 x 10^ fs^ and a — 1 radian, re- 
spectively. Black crosses in (b) represent the observed signal 
for a random (instead of quadratic) distribution of phases 
in the pulse train. In (c) and (d), efficiency of the population 
transfer is shown as a function of the continuous and piecewise 
chirp, respectively. In both cases, the energy of the excitation 
field corresponds to the first minimum of the respective Rabi 
oscillation. Experimental signal on all 4 plots is normalized 
to the maximum of Rabi oscillations in plot (a). 



transform-limited pulse. This and all other experimen- 
tal data points were normalized to the first maximum of 
this curve at around 0.1 /iJ. In plot (b), the similar os- 
cillatory dependence corresponds to the train of 7 pulses 
with Gaussian envelope of amplitudes and zero piecewise 
chirp, namely Ak = exp(— /c^ In 2/4) and <l>/e = where 
—3 < < 3. The temporal separation between the pulses 
in the train was 400 fs. As expected, the integrated pulse 
area of the train is larger than the corresponding area of 
a single pulse with the same total energy. The difference 
stems from the different scaling of the pulse area and 
pulse energy with the field amplitude (first and second 
power, respectively), resulting in the smaller period of 
Rabi oscillations for the piecewise excitation. 

When either continuous or piecewise chirp is applied 
to a single excitation pulse or a series of pulses, respec- 
tively, the amount of excited Rb atoms ceases to oscillate 
with the pulse energy. In both cases it saturates at the 
maximum value and becomes insensitive to the excitation 
field strength in the typical AP way (red circles in Fig. 
|4]^a,b)). We attribute the decay of both the amplitude 
of the Rabi oscillations and the saturated signal to the 
weak pre-pulse generated by our laser system 3 ns prior 



to the main pulse. To demonstrate the significance of the 
quadratic phase in the piecewise excitation, in Fig. |4]^b) 
we also plot the observed signal for a pulse train with 
random (black crosses). Even though the mutual co- 
herence between the pulses is preserved regardless of the 
random phase jumps, such pulse sequences do not result 
in the AP-like saturation of the population of the excited 
state. 

The stability of the adiabatic passage can be further 
seen in the dependence of the transfer efficiency on the 
magnitude of the chirp. In both the continuous and piece- 
wise scheme, we set the unchirped pulse area to 27r, which 
corresponded to the energy of 0.31 and 0.07 /iJ, respec- 
tively. The result of scanning the conventional chirp from 
— 10^ to +10^ fs^ is similar to scanning the piecewise 
chirp a between -1.5 and +1.5 radians, attesting to the 
similar mechanisms of the two processes. Note the decay 
of the excited state population at \a\ > 1.5 rad in plot 
(d), caused by the breakdown of the piecewise adiabatic- 
ity due to the increasingly high phase increments from 
pulse to pulse in the pulse train. The results of numeri- 
cal simulations are shown along with the corresponding 
experimental results as solid lines in Fig. [4j The only 
single fitting parameter for all calculations was the area 
of the excitation beam. FROG traces were used to define 
the temporal profile of the pulses. 

In contrast to adiabatic transfer with continuously 
chirped pulses, PAP can be easily generalized to a wide 
class of cases. For example, one can consider replacing 
states |1) and |2) by wave packets composed of many indi- 
vidual eigenstates. Piecewise excitation of such systems 
could be advantageous because, by their very nature, 
wave packets often evolve periodically in time. Viewed 
from the spectral perspective, this extension of PAP is 
related to the "Coherently Controlled Adiabatic Pas- 
sage" schemes, proposed earlier for controlling popula- 
tion transfer between wave packets [16j. We have simu- 
lated numerically a number of situations in which states 
|1) and 1 2) correspond to an isolated vibrational level 
and a wave packet composed of several vibrational levels, 
respectively. The population transfer is executed with 
a train of femtosecond pulses separated in time by the 
wave packet's vibrational period. Quadratically vary- 
ing the carrier phase of the pulses enables robust and 
complete population transfer into the excited state wave 
packet, whose shape can be controlled by spectrally shap- 
ing the individual pulses in the train. The availability of 
such robust and selective population transfer in multi- 
level systems positions PAP as a powerful tool in con- 
trolling molecular dynamics. 
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